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Abstract. X-Ray observations of groups and clusters of galaxies show that the Intra-Cluster Medium (ICM) in 
their cores is hotter than expected from cosmological numerical simulations of cluster formation which include star 
formation, radiative cooling and SN feedback. We investigate the effect of the injection of supersonic AGN jets into 
the ICM using axisymmetric hydrodynamical numerical simulations. A simple model for the ICM, describing the 
radial properties of gas and the gravitational potential in cosmological N-Body+SPH simulations of one cluster 
and three groups of galaxies at redshift z = 0, is obtained and used to set the environment in which the jets 
are injected. We varied the kinetic power of the jet and the emission-weighted X-Ray temperature of the ICM. 
The jets transfer their energy to the ICM mainly by the effects of their terminal shocks. A high fraction of the 
injected energy can be deposited through irreversible processes in the cluster gas, up to 75% in our simulations. 
We show how one single, powerful jet can reconcile the predicted X-Ray properties of small groups, e.g. the 
Lx — Tx relation, with observations. We argue that the interaction between AGN jets and galaxy groups and 
cluster atmospheres is a viable feedback mechanism. 
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1. Introduction 

Clusters and groups of galaxies are observed in the X-ray 
band mainly due to the thermal bremsstrahlung radiation 
emitted by the hot gas, the Intra Cluster Medium (ICM). 
A wealth of recent observations, coming in particular 
from the X-ray telescopes CHANDRA and XMM-Newton 
(Paerels & Kahan 2003), have evidenced a complex ICM 
structure, which is characterized not only by cooling cores 
(Peterson et al. 2003, Kaastra et al. 2004), but also by cold 
fronts, bow-shocks, "bullet-like" structures (Vikhlinin et 
al. 2001, Markevitch et al. 2002, Kempner et al. 2003). 
While all of these features can be understood in terms of 
the currently favored hierarchical model of structure for- 
mation, there is an ongoing debate concerning the origin 
and the details of the energy budget of the ICM. 

In fact, the most obvious way to explain the heating of 
the ICM to the temperatures observed in groups and clus- 
ters of galaxies, ranging from « 0.5 to more than 10 keV, 
is to suppose a pure gravitational heating. Such a heat- 
ing is a natural outcome of the merging processes which 
form groups and clusters of galaxies in the hierarchical 
model picture. While the bigger and bigger Dark Matter 
halos merge, their gas is shocked and heats up, to form the 
groups and clusters atmospheres which are observed. The 
power spectrum of the initial density perturbations grow 



hierarchically to form groups and clusters and does not 
possess a characteristic length or mass scale, in the rele- 
vant range. The gravity itself has no characteristic scale. 
Thus, with the additional assumption that gas is in hy- 
drostatic equilibrium with the cluster potential, it can be 
shown that this scenario predicts simple self-similar scal- 
ing relations among the various observable X-ray proper- 
ties of groups and clusters. In particular, the relation be- 
tween X-ray luminosity and the temperature is Lx oc T^, 
and the entropy S = T jn'J'^ (where is the electron 
number density of the gas) scale as S" oc T when S is 
evaluated at fixed overdensities for different clusters. 

Nevertheless, the observational situation appears to 
be quite different. The luminosity-temperature relation 
is steeper than this simple picture predicts, Lx cc T" 
with a ~ 2.5 — 3 for clusters with temperature T > 2 
keV (White, Jones & Forman 1997; Markevitch 1998; 
Arnaud & Evrard 1999; Ettori, De Grandi & Molendi 
2002). There are indications of a steeper slope for groups 
with T < IkeV (Ponman et al. 1996, Helsdon & Ponman 
2000, Sanderson et al. 2003). The gas entropy in clusters 
is found to be in excess with respect to the model predic- 
tions (Ponman, Cannon & Navarro 1999; Lloyd-Davies, 
Ponman & Cannon 2000; Finoguenov et al. 2003), and its 
dependence on the temperature approximately follows the 
law S oc T^^^ (Ponman, Sanderson & Finoguenov 2003). 
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The discrepancy between observations and theory sug- 
gests that the ICM is subject to a heating which is larger 
than what the simple self-similar model predicts. In gen- 
eral, there is a consensus on the fact that the observations 
can be matched using models where a given amount of 
energy per particle is added to that coming from gravi- 
tational processes, independently on the precise origin of 
this energy. Numerical studies suggest that a heating en- 
ergy of about 0.5 — 1 kcV per particle, or cquivalcntly 
an entropy floor of « 50 — 100 keV cm^, is enough to 
satisfy the observational constraints (Bialek et al. 2001, 
Brighcnti & Mathews 2001, Borgani ct al. 2001, 2002). 
Semi-analytical models give an even higher estimate of 
the needed heating, up to ~ 400 keV cm^ (Babul et al. 
2002, McCarthy, Babul & Balogh 2002). The ICM can be 
heated by two natural processes, which are present when 
the baryons physics is considered. The first is the energy 
feedback from supernovac explosions, which directly heats 
cold gas and increases its entropy. The second is the radia- 
tive cooling of the gas itself, which subtracts low-entropy 
cold gas from the hot phase and turn it into stars, thus 
increasing the average temperature and entropy of the re- 
maining gaseous medium. The efficiency of the two pro- 
cesses in heating the ICM have been extensively studied 
so far, both with semi-analytical techniques (Cavaliere, 
Menci & Tozzi 1998; Tozzi & Norman 2001; Menci & 
Cavaliere 2000; Bower et al. 2001; Voit & Bryan 2001; 
Volt et al. 2002; Wu & Xue 2002) and with numerical 
simulations (Evrard & Henry 1991; Kaiser 1991; Bower 
1997; Balog et al. 1999; Muanwong et al. 2002; Dave et al. 
2002; Babul ct al. 2002). 

Recent numerical studies also considered simultane- 
ously the effect of radiative cooling, star formation and 
feedback from SN (Suginohara & Ostriker 1998; Lewis et 
al. 2000; Yoshida ct al. 2002; Loken et al. 2002; Tornatore 
et al. 2003; Borgani et al. 2004). These studies partially 
disagree on the resulting agreement with observations 
which can be obtained, ranging from a more (Loken et 
al. 2002) to a less (Borgani ct al. 2004) optimistic inter- 
pretation of the simulations results. All of them, however, 
showed how the fraction /, of gas which is converted into 
a "cold" stellar phase results to be too high, in excess with 
respect to the value 10% which is obtained from mea- 
surements of the local luminosity density of stars (Balog 
et al. 2001; Lin, Mohr & Stanford 2003). The overcoohng 
problem which naturally descend from these results still 
calls for the need of further energy contributions to the 
ICM to be solved. 

A possible source for this energy is given by AGNs 
activity. The available energy budget is in principle more 
than adequate to heat the ICM to the desired temperature 
(Valageas & Silk 1999, Cavaliere, Lapi & Menci 2002). In 
fact, the accretion of gas onto supermassive BHs in galac- 
tic cores gives outputs of order 2 • 1O^^(M/1O^M0) erg, 
where M is the accreted mass and a standard mass-energy 
conversion efficiency of 10^^ is used (Wu, Fabian & Nulsen 
2000; Bower et al. 2001). Observations with the ROSAT 
satellite suggested an interaction of the central AGN ac- 



tivity with the ICM in the Perseus cluster (Boringer et al. 
1993), possibly under the form of "bubbles" (McNamara, 
O'Connel & Sarazin 1996). CHANDRA also gave obser- 
vational hints of possible coupling between ICM and cen- 
tral radio-sources in clusters: observations of the Perseus 
cluster by Fabian et al. (2000) showed that the cluster gas 
is displaced by the radio lobes of the source 3C84 which 
is inflating cavities or "bubbles" in the X-ray emitting 
gas; observations of ripples and sound waves propagat- 
ing through the ICM (Fabian et al. 2003a) suggest that 
the cavities are gently expanding and interacting with the 
surrounding medium. A similar scenario hes been depicted 
in the case of M87 by Forman et al. (2003) who observed 
weak shock fronts associated with the AGN outburst prop- 
agating through the surrounding gas. In the case of the 
Hydra A cluster, hosting the powerful FRI radio source 
3C218, it has been argued that, even if there are no signs 
of shock-heated gas around the radio lobes (McNamara 
et al. 2000), shock heating is needed in the cluster core 
to balance cooling flows (David et al. 2001). Recent deep 
CHANDRA observations of the same source (Nulsen et al. 
2004) have revealed features in the X-Ray bightness pro- 
file which can be associated with a shock front driven by 
the expanding radio lobes requiring a total energy around 
]^q61 g]-g drive it. Observations of Cygnus A (Smith et al. 

2002) have revealed signs of stronger interaction between 
the radio source and the ICM since the shells of enhanced 
X-ray emission around the cavity are hotter than the sur- 
rounding medium perhaps as a result of heating by a bow- 
shock driven by the expanding cavity. X-ray emission from 
shocked gas has been detected in the observations of the 
radio galaxy PKS 1138-262 at z=2.156 (Carilh et al. 2002) 
for which a jet kinetic power 5 x 10^® erg s^^ injected for 
2 X 10*" years has been estimated. The X-ray emission ob- 
served around the FRII radio source 3C294 at z=1.786 
(Fabian et al. 2003), if interpreted as thermal radiation 
from shocked material, gives an estimates of the injected 
power greater than 10^^ erg s~^. 

On the other hand, we still miss a coherent and physi- 
cally well based model of the coupling of this energy with 
the ICM. Different studies of coupling between the AGN 
power and the ICM have been done, which can be roughly 
divided in two main classes: the study of the buoyant be- 
havior of hot gas bubbles in the ICM, almost indepen- 
dently of how the bubbles were produced (Quills, Bower 
& Balogh 2001; Churazov et al. 2001, 2002; Briiggen et al. 
2002; Briiggen & Kaiser 2002; De Young 2002; McCarthy 
et al. 2003; Mathews et al. 2003; Hoeft, Briiggen & Yepes 

2003) , and the direct study of the effect of injection of me- 
chanical energy, via an AGN jet, in the ICM (Cavaliere, 
Lapi & Menci 2002; Reynolds, Heinz & Begelman 2002; 
Nath & Roychowdhury 2002). In the first view, the en- 
ergy supply to the cluster is gentle and buoyancy-driven; 
in the second view, it is a momentum-driven process, char- 
acterized by the formation and dissipation of shock fronts. 
Other ideas have been proposed as well, for instance the 
viscous dissipation of sound waves (Pringle 1989) which 
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can arise from the activity of a central AGN (Ruszkowski, 
Briiggen & Begelman 2003). 

Semi-analytical studies of the above processes are usu- 
ally based on simple prescription for the exchange of en- 
ergy between the central AGN and the ICM, coupled with 
basic models for describing the dynamical evolution of 
bubbles, cocoons or jets (Gavaliere, Lapi & Menci 2002; 
Nath & Roychowdhury 2002; Churazov et al. 2002; De 
Young 2002; McCarthy et al. 2003; Mathews et al. 2003) 
and, at best, a modeling of the cluster atmosphere (Wu, 
Fabian & Nulsen 2000) , usually consisting either in a clas- 
sical /3-model (Gavaliere & Fusco-Femiano 1976) or in a 
modified model of the same kind, in which the hot clus- 
ter gas in in hydrostatic equilibrium with the gravita- 
tional potential of the galaxy cluster. In the last case, the 
mass density profile of the cluster is usually taken to be a 
Navarro, Frenk & White (NFW) analytical profile, which 
have been shown to correctly describe, on average, the 
shape of such clusters in cosmological N-body simulations 
(Navarro, Frenk & White 1996, 1997). The gravitational 
potential is then derived consequently. 

Hydrodynamical numerical simulations of the IGM 
heating by AGN also resort to these class of models for 
the cluster atmosphere, sometimes slightly modified for 
closely following an observed distribution of the mass and 
electron density temperature (Nulsen & Bohringher 1995). 
Simulations of buoyant bubbles are usually performed in 
three dimensions, at the cost of a lower resolution (Quilis, 
Bower & Balogh 2001; Briiggen et al. 2002) or in two di- 
mension assuming spherical symmetry thus reaching good 
resolution only in the central part of the computational 
domain (Churazov et al. 2001). 

Simulations of jets in the ICM have been performed 
in 2D with spherical symmetry by Reynolds, Heinz & 
Begelman (2002). Their jet was injected in a cluster at- 
mosphere for w 5 ■ 10'' yr. The atmosphere was taken to 
be a simple /3-model, and a suitable static gravitational 
potential was added to guarantee the hydrostatic equi- 
librium. No attempt has been done there to model the 
differences which can arise in the ICM as a function of 
the group/cluster mass, nor to have a mass profile simi- 
lar to the one coming from cosmological simulations. The 
atmosphere was then evolved for as long as 10^ yr. The 
main result of this work is that, when an equilibrium is 
achieved again by the cluster atmosphere, its core specific 
entropy has been enhanced by ~ 20 per cent and a large 
fraction of the injected energy, as high as « 50 per cent, is 
thermalized onto the ICM. A three dimensional study of 
the interaction between a low power jet and the ICM has 
been performed by Omma et al. (2004) having a resolution 
around w 0.6 kpc with 1024 cells per side. 

A numerical study of viscous dissipation of waves had 
still higher length resolution (« 0.1 kpc with 2048 cells 
per side, Ruszkowski, Briiggen & Begelman 2003). They 
showed how the viscous dissipation of the energy of sound- 
waves produced by a single AGN duty-cycle is insufhcient 
to balance radiative cooling of the gas; but intermittent 



activity of the central source is able to offset cooling by 
this process. 

Finally, a cosmological self-consistent study of the evo- 
lution of the radio plasma during a merger of cluster of 
galaxies has been performed by Hoeft, Briiggen & Yepes 
(2003). They used the Tree-I-SPH N-body code GADGET 
(Springel, Yoshida & White 2001) and achieved a high 
mass and force resolution (the mass of a gas particles was 
Mgas = 2.6- lO^/i~^M0 and the Plummer-equivalent soft- 
ening length was e = 2h~^ kpc). However, their focus was 
on the evolution of the fossil radio-relics and not on the 
energy exchanges between the AGN which produced such 
relics and the ICM. 

In the present work, we study the energy exchange 
between AGN-powered jets and the surrounding ICM of 
group/cluster of galaxies, with 2D hydrodynamical nu- 
merical simulations. We extracted the groups and clus- 
ters atmospheres from cosmological numerical simulations 
(Tornatore et al. 2003), at the redshift z = 0. We build 
a self-similar model for such atmospheres based on the 
NFW radial density profile for the mass distribution and 
the hydrostatic equilibrium with the consequent potential 
for the gas. We note that our model is not an average de- 
scription of a typical group or cluster of galaxy, but it has 
been built to represent the output of a self-consistent nu- 
merical cosmological simulation. The jet is injected in the 
ICM for a time t ^ 2 - 10^ yr, then switched off. We varied 
the temperature scale of our model and the kinetic energy 
of the jet. We derive the energy balance of the ICM, the 
entropy evolution, and, for the first time using this ap- 
proach, the resulting Lx — Tx relation, after the heating 
of the ICM has occurred. 

The plan of the paper is the following: in Section 12 
we describe the baseline model. In Section |31 the numeri- 
cal simulations are presented. Section 01 is devoted to the 
calculation of the energy balance and Section [3 to the 
entropy evolution. In SectionElwe calculate a new hydro- 
static equilibrium for the ICM and derive the Lx — Tx, 
Lx — and S ~Tx relations. In Section|Slwe draw our 
conclusions and discuss future perspectives of this work. 

2. The baseline model 

A self-consistent, numerical study of the interplay between 
AGN jets and clusters atmospheres would require a huge 
dynamical range, to resolve simultaneously the kpc scale 
of the jets and the « 100 Mpc scale required for following 
the clusters formation and evolution. Such a dynamical 
range is still outside the current computing power, even 
for supercomputers. Moreover, the details of the formation 
and evolution of the AGN itself are not completely un- 
derstood, thus they cannot be parametrized easily in self- 
consistent, cosmological numerical simulations. To achieve 
high spatial resolution, we perform 2D hydrodynamical 
simulations in cylindrical symmetry. Therefore, we can- 
not directly use the N-Body-|-SPH simulated cluster at- 
mospheres as initial conditions: being 3D, their resolu- 
tion is significantly lower than the one we are going to 



4 



Zanni et al.: Heating groups and clusters of galaxies: the role of AGN jets 



reach. We thus build a simple model for the cluster atmo- 
spheres, to reproduce the spherically averaged properties 
of the cluster gas in cosmological Tree-|-SPH simulations. 
We consider adiabatic simulations performed by Tornatore 
et al. (2003, T03 hereinafter) of the evolution of four clus- 
ters having virial mass Mvir = 2.35, 2.52, 5.98, 39.4 x 
10^"^ Mq. We take the clusters at z = in order to di- 
rectly estimate the effect of jets energy injection on the 
X-ray properties of galaxy clusters and to compare these 
properties with observations. The parameters of our model 
are constrained using the radial profiles of these simulated 
objects. An additional requirement is that the model does 
scale with cluster mass. Our initial conditions for the en- 
vironment in which the jets propagate are then set, in our 
hydrodynamical simulation, using such model. As in T03, 
we use a cosmology with fi^ = 0.3, fi^ = 0.7, Hubble 
constant i?o = 100 h km s^^ Mpc with h — 0.7, and 
^^bar = 0.019 which gives a cosmic baryon fraction 
of /bar — 0.13. The baryonic gas is taken with a primor- 
dial composition (hydrogen mass- fraction X = 0.76 and 
Y = 0.24 for helium, mean molecular weight fj, = 0.6). We 
define here the virial radius i?vir as the radius of a spher- 
ical volume inside which the mean density is Ac(z) times 
the critical density (Afvir(-2) — 47ri?^ij,pcrit Ac/3), where 
Ac{z) is the solution of the spherical collapse of a top-hat 
perturbation in an expanding universe (Bryan & Norman, 
1998). At z = 0, for our chosen cosmology is Ac = 101. 

Firstly we model the profiles of the dark matter den- 
sity taken from the T03 simulations with a NFW profile 
(Navarro, Frenk & White 1996): 



Pdii 



PdmO 



x{l 



(1) 



with X = r/vs where rg is a characteristic radius. The 
dark matter central density pdmO is determined by requir- 
ing that the dark matter mass contained inside the virial 
radius is equal to (1 — /bar)Afvir- The only free parameter 
of the dark matter distribution is therefore the concentra- 
tion parameter c = Rviv/rg defined by the ratio between 
the virial radius of the cluster i?vir and the characteris- 
tic radius r^: from the T03 simulations, we estimate an 
average value of c = 9.5 for our four objects. We use 
this value of the concentration parameter for modeling 
the DM profiles of every cluster, since we want to obtain 
a model whose properties do scale with the virial radius 
of the objects (via the value of the DM central density). 
This assumption is not completely correct, since it is well 
known (Navarro, Frenk & White 1996, 1997; Bullock et al. 
2001) that the concentration of a dark matter halo is de- 
termined by its mass and its formation history; at a given 
mass, early formation generally leads to a denser core and 
thus to a higher concentration. However, taking the aver- 
age value for c from the Tree-|-SPH simulations at z = 
we have a maximum error of 10% on the individual values 
of the concentration parameter. 

In our calculations, the distribution of the dark matter 
density determines the shape of the gravitational poten- 
tial well. This potential will be considered steady in time: 



this approximation is reasonable, since the cosmological 
timescale for the evolution of a DM halo is much longer 
than the dynamical timescale of an AGN jet. The gravi- 
tational potential cj) generated by the density distribution 
of Eq. (P) is 



G(l - /bar)Afvir (1 + x) 



Rv 



log(l+e) 



1 

1+c 



(2) 



The gas distribution is determined imposing hydrostatic 
equilibrium with the potential of Eq. ^ and neglecting 
the self gravity of the gas. Assuming a polytropic equation 
of state for the gas in initial equilibrium {P — Pa{p/ po)^) 
we obtain 



1 



log(l + x) 



1 



i/(r-i) 

n (3) 



The value taken for the polytropic index is F = 1.1, in 
order to reproduce the decrease of temperature observed 
in the T03 simulated clusters and groups at z = 0. The 
coefficient (3 given by the expression 



/3 = 



G(l - /bar)Afv 



7 



Rv 



log(l-l-c) ]_ 

c 1+c 



(4) 



is taken constant at all the physical scales in order to 
have a self-similar density profile and is equal to f3 — 8, 
which is a good approximation for the T03 simulations. In 
Eq. 10} the index 7 is the ratio of the specific heats and 
Cso = -y/ jkTo I pm^i is the central sound speed. Finally the 
gas central density po is determined by imposing that the 
gas mass inside the virial radius is equal to /bar Mvir. 

With all the above assumptions we obtain the following 
scaling relations for our baseline model: 



nco = 0.033 cm 



KTo = 7.43 



i^Tow = 6 



2/3 



keV 



IOI^Mq 



2/3 



keV 



KTm^ = 4.5 



n = 272 



1O15M0 



2/3 



keV 



(5) 



1/3 



kpc 



/ KT 



Lb,. = 4.34 X 10« erg s'' 



V keV 



<5'o. 



1 -R21 



2/3 



36.3 



KT 
keV 



keV cm 



for, respectively, the central electron density Uco, the cen- 
tral temperature KTq, the emission weighted temperature 
KTcw, the mass weighted temperature KT^^, the charac- 
teristic radius Vg, the bremsstrahlung luminosity Lbr and 
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the entropy S calculated at 0.1 i?200 as in Ponman et 
al. (2003). The emission weighted temperature, assuming 
bremsstrahlung emissivity, and the mass weighted tem- 
perature are defined as 

r-Rvi, 



Jo 



C"pTdV 
C" Pdv 



Anyway, in the following, we will also include metal-line 
cooling to compute the X-Ray luminosities. Our assump- 
tion of a constant concentration parameter for all clusters 
does reflect in a constant central electron density. Since 
our baseline model is self-similar at all scales, as it is 
shown by Eqs. lO, all the characteristic parameters de- 
scribing the model (central density, central temperature, 
virial mass, etc.) scale with only one free parameter. Such 
parameter can be chosen as the mass Mvir or equivalently 
the virial radius i?vir, the temperature Tew, etc., of the 
cluster or group. 

Having derived an analytical expressions for the dis- 
tribution of all the physically interesting quantities and 
determined the parameters from the T03 Tree-|-SPH cos- 
mological simulations, we can finally test how well the 
model agrees with the radial profiles, i.e. how accurately 
the T03 data are reproduced. 

In Figs. Hand 121 we show the comparison between the 
radial dependencies of some quantities as derived from 
the T03 simulations and as modeled by us for two clus- 
ters of galaxies with respectively a virial mass of Mvir = 
39.4 X 10^3 j^^.^ ^ 2.52 x 10^^ Mq. The agree- 

ment of the model with the simulation data appears to 
be good for both objects, as far as p, pdm, ^Tot, P, ans 
iS* are concerned. There are slight discrepancies only in 
the behavior of T, especially for the most massive clus- 
ter of Fig. ^ where the central drop is not captured by 
the model (see however T03 for a discussion of tempera- 
ture profiles in SPH simulations). As a consequence, the 
analytical model also reproduces the behavior of the X- 
ray properties shown by T03 numerical simulations. We 
therefore conclude that the initial conditions for our hy- 
drodynamical simulations can safely been set using the 
model presented in this Section. 

3. The simulations 

In order to investigate the interaction of an extragalac- 
tic jet with the intracluster gas, we performed two di- 
mensional hydrodynamic numerical simulations of a cylin- 
drical (axisymmetric) jet propagating in a gravitationally 
stratified but not self-gravitating atmosphere: 

dp 

dt 

dv 

— + (t) • S/)v + Vp/p = V0 



dp 
di 



dt 

+ {v ■ V)p - 7 



(6) 



P 



dp 
dt 



+ {v ■ V)p 







where the fluid variables p, p and v are pressure, den- 
sity and velocity respectively; 7 = 5/3 is the ratio of the 
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Fig. 1. Comparison between our baseline model {dashed 
line) and and the results of the adiabatic simulation by 
Tornatore et al. (2003) {solid line) of a cluster with a mass 
Mvir = 39.4 X 10^3 Mq, taken at z = 0. We show the radial 
dependencies for the following quantities: ( Upper panels) 
gas density p, dark matter density pdm: total mass Mtot; 
{Lower panels) gas pressure P, gas temperature T and gas 
entropy S = T/ nj'^ 
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Fig. 2. Same comparison as Fig. ^ for a group with a 
mass Mvir = 2.52 x 10^^ Mr... 



speciflc heats. Radiative losses are neglected. The cooling 
times in the core of the unperturbed clusters are longer 
than the time of activity of our jets (see below). The cool- 
ing times in the centers are of the same order of magnitude 
of the elapsed time at the end of the simulations. On the 
other hand, the material near the center of the clusters is 
subject to signiflcant expansion due to the energy injection 
itself and therefore we expect the importance of cooling to 
be reduced. The system of Eqs. (O has been solved numer- 
ically employing a PPM (Piecewise Parabolic Method) hy- 
drocode (Woodward & Colella 1984). As discussed in Sect. 
Elthe density profile of the undisturbed ambient medium is 
given by Eq. Q while the gravitational potential needed 
to take the gas distribution in equilibrium is given by Eq. 
. Simulations have been performed to test the numerical 
stability of the initial gravitational equilibrium revealing 
that small subsonic speeds begin to appear in the central 
core of the initial profile on timescales much longer than 
the energy injection time thus ruling out the possibility 
of spurious heating effects. A (cylindrical) jet is injected 
from the bottom boundary of the integration domain, in 
pressure balance with the ambient. The axis of the jet is 
along the left boundary of the domain (r = 0), where we 
have imposed symmetric boundary conditions for p, p, 
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and antisymmetric conditions for Vr- Reflective boundary 
conditions are also imposed on the boundary of injection 
of the jet (z = 0) outside its radius in order to reproduce 
a bipolar flow and to avoid spurious inflow effects. On the 
remaining outer boundaries we extrapolated the equiUb- 
rium profiles in order to keep the initial atmosphere in 
equilibrium on the outer parts of the computational do- 
main. 

Measuring lengths in units of the characteristic radius 
Ts, velocities in units of the central adiabatic sound speed 
CsO in the undisturbed external medium and the density in 
units of the ambient central density po, our main param- 
eters are the Mach number M = v-Jcso and the density 
ratio v = Pj/ Po- In the following discussion we will take 
the emission weighted temperature T^^i as free parameter 
for describing the unperturbed medium. Consistently the 
unit of time is 



<o = — = 1.9 X 10" years 

CsO 

Defining the units for the jet kinetic power as 



(7) 



= -zPor^ c; 



2„3 

i "^so 



where rj is the jet radius, the kinetic power Lj can be 
expressed as 

Lj = APvL], . (8) 
Assuming a jet radius r-^ = A kpc ik is given by 



3/2 

Lk = 3 X 10^" ( ) erg s~i 



-,42 / ^ew \ 

IkeVy 



(9) 



The cylindrical symmetry adopted in our simulations, con- 
trary to spherical symmetry, allows us to have a high reso- 
lution from the jet scale up to the cluster scale. We inject 
a jet having 0° opening angle. The assumed jet radius, 
dictated mostly by the need of an adequate resolution (at 
least 10 grid points) on the jet radius scale, is consistent 
with the size of jets with an opening angle around 9° — 15° 
(Feretti et al. 1999, Lara et al. 2004) at a distance ~ 15-25 
kpc from the parent core. Moreover we must consider that 
the lifetime of the jets is at least one order of magnitude 
shorter than the total time followed in the simulations 
and after the jet has ceased its activity, as we will show in 
detail in the next Section, it loses its coUimation and com- 
pletely disappears: from that time onward the assumed jet 
radius is no more important and the whole dynamics of 
the system is determined, as we will show later, by the 
total energy which has been injected. In order to follow 
separately the evolution of the injected material and the 
ambient medium we define a set of four passive tracers 
following the evolution equation 



dt 



+ [v ■ V)7I 



The tracers are initialized as follows: 71 = 1 for r/r^ G 
[0,1], T2 = 1 for r/rs G ]1,3], % = 1 for r/r, G ]3,6], 
74 = 1 for r/fs G ]6,9.5] and equal to otherwise. 



We perform a first set of six simulations for three 
different X-ray temperatures (Tew — 0.5, 1, 2 keV), 
which correspond in our baseline model to virial masses 
Mvir = 2.3 X 10^3^ 6.69 x lO^^, 1.89 x lO" Mq respec- 
tively, and with two different jet kinetic powers (Lj = 10**^, 
5 X 10**^ erg s""'^). The masses of the clusters used as ini- 
tial conditions for our simulations are different from the 
virial masses of T03; we used T03 simulations only for 
building our baseline model. In order to evaluate the ef- 
fects on the large scale structures of an identical source 
of energy, the jets of equal power injected in the three 
groups/clusters have the same characteristic parameters 
(i.e. density, speed and radius) and the two different pow- 
ers are determined by their jet speed. Since in our model 
the central density of the ambient medium is a constant at 
all scales we choose then a constant v = 0.001 for all these 
simulations. Our simulated cases are then characterized by 
their respective Mach number. In all of our simulations we 
choose a square computational domain with a resolution 
of 250 points for 100 kpc: the size of the computational 
domain ^domain is determined by the minimum between 
the virial radius of the structure and the maximum length 
scale that we decided to simulate (819.2 kpc or equiva- 
lently 2048 points). The details about the computational 
domain can be found in Table ^ In this set of simulations 
the jet is injected into the computational domain for a 
time iiife = 0.1 to — 2 x 10^ years and is then switched 
off imposing reflecting boundary conditions over the whole 
boundary of injection (z = 0). In this way the total energy 
injected into the computational domain is of the order of 
i?j ~ 1.2 x 10^^ or i?j ~ 6 X 10^^ erg depending on the 
jet power considered and taking into account the energy 
injected by both the jet and the counterjet. This energy 
has to compete with the huge thermal energy of the group 
or cluster i?th which for our baseline model is given by 



3.1 X lO^^i ( — 
\ kcV 



5/2 



erg 



(10) 



After the jets have ceased their activity at i/^o = 0.1 the 
simulations are followed until the expanding bow shock 
determined by the energy injection covers at least 80% of 
the computational domain. In Table^we show the param- 
eter M for the simulations that we have performed giving 
also the values of the ambient temperature Tqw, the en- 
ergy injected into the domain £j in units of the thermal 
energy Eth of the cluster, the size of the domain ^domain, 
the ratio between the size of the domain and the virial 
radius -Rvir, the number of grid points ridomain used for 
every of the two dimensions and the final simulated time 
given in units of to- 

In order to understand how much the results are sensi- 
ble to the jet parameters considered and to the resolution 
used, we also add three test simulations: using as initial 
conditions an unperturbed atmosphere with a mean tem- 
perature Tew — 0.5 keV, we first repeat the case with 
Tj = 10*6 erg s-i and M = 211 with hah the spatial 
resolution. Then we consider a jet with a kinetic power 



Zanni et al.: Heating groups and clusters of galaxies: the role of AGN jets 



7 



Table 1. Parameters of the six simulations performed: Mach number M, emission weighted temperature of the ambient 
Tew, jet power Lj, injected energy Ej/Eth, size of the domain /domain, size of the domain in unit of the virial radius 
'domain /^vir> number of grid points ^domain and final time of the simulations rfi„ — t/tQ. 



M 


Tow (keV) 


Lj (erg s ^) 


Ej/Eth 


^domain (kpc) 


^domain / Rvir 


^domain 


Tfln 


105 


2 


10* 


0.068 


819.2 


0.55 


2048 


3.8 


150 


1 


10* 


0.39 


819.2 


0.77 


2048 


3.9 


211 


0.5 


10* 


2.2 


744.5 


1 


1860 


2 


180 


2 


5 X 10* 


0.34 


819.2 


0.55 


2048 


2 


255 


1 


5 X 10* 


1.9 


819.2 


0.77 


2048 


1.95 


361 


0.5 


5 X 10* 


10.9 


744.5 


1 


1860 


0.85 



Table 2. Parameters of the three test simulations per- 
formed: Mach number M, density ratio i' = Pj/po, ki- 
netic jet power Lj , active phase iufc, number of grid points 
'T'domain and final time of the simulations Tfin = t/tQ. All 
these simulations have as initial condition an atmosphere 
with Tew = 0.5 keV. 



M 




Lj (erg s ^) 


tiifc (10' yr) 


^domain 


Tfln 


211 


0.001 


10* 


2 


930 


3 


361 


0.0002 


10* 


2 


1860 


3 


211 


0.0005 


5 X 10* 


4 


1860 


2 



U = 10 



46 gj.g g 1 injected for t\ 



ife = O.lto — 2 X 10^ years 
as before but with a different Mach number and density 
ratio. Finally we simulate a jet with a lower kinetic power 
(lower density) injected for a longer time, taking the total 
injected energy equal to Ej ~ 1.2 x 10^^ erg as before. In 
Table[51wc show the parameters of these three simulations. 



4. The energy balance 

Firstly, we need to understand how the extragalactic jets 
interact with the surrounding medium and how the energy 
injected by the jet couples with it. Theoretical models (see 
Begelman and Cioffi 1989) and numerical simulations (see 
Norman et al. 1982, Massaglia et al. 1996) have analyzed 
in detail the structure of the interaction between an un- 
derdense supersonic jet and the denser ambient: the jet is 
slowed down by a strong terminal shock (Mach disk) where 
its material is thermalized and inflates an expanding co- 
coon that compresses and heats the ambient medium. The 
compressed ambient medium forms a shell surrounding the 
cocoon which is separated from the shell itself by a con- 
tact discontinuity. During this phase the cocoon and the 
shell are usually strongly overpressured with respect to the 
outer pressure (i.e. the initial ambient energy enclosed in- 
side the shell's outer radius is negligible compared to the 
injected energy): the shell of compressed ambient material 
expands therefore with a supersonic speed driving a strong 
shock into the surrounding medium where the expansion 
work done by the cocoon of jet material is irreversibly 
dissipated. In this regime the temporal behavior of the 
expanding shock can be obtained qualitatively assuming 



spherical symmetry (that is a good approximation at least 
for light jets; see Zanni et al. 2003 and Krause 2003) and 
the following behaviors for the injected energy Ej and the 
density p stratification: 



Ej - Ea 



Pa 



(11) 



Since there are no characteristic lengths defined, it is pos- 
sible to derive the temporal behavior of the expanding 
Shockwave only from dimensional arguments (see Sedov 
1959): we can write the shock Rg radius as 



Rs oc 



Eorl 
Pott 



d + 2 



(12) 



Assuming a polytropic relation for the ambient pressure 
P (X p'^ we can derive the following behaviors for the shock 
speed Vs and Mach number Ms'. 



and 



kd(l-7) + 2d-2-yk-6 

M, OC t 



(13) 

It is easy to show that the Mach number decreases with 
time if the rate at which the energy is injected (d) is lower 
than the rate at which the background energy enclosed 
in the shock radius (oc J^^" Pr^ dr) grows: this means 
that as the shock expands the ambient energy enclosed 
inside the shock volume becomes comparable to the in- 
jected energy and so the cocoon slowly reaches a pressure 
equilibrium with the surrounding atmosphere expanding 
almost at the local sound speed. The contact discontinu- 
ity, which in the overpressured phase expands at the same 
rate of the outer shell radius (see Zanni et al. 2003), is sub- 
ject to Rayleigh- Taylor instabilities if its deceleration does 
not exceed the local gravity acceleration and to Kelvin- 
Helmholtz instabilities: these effects give rise to mixing 
between the jet material which forms the cocoon and the 
shell which surrounds it. Obviously, the entrainment with 
the high entropy material of the cocoon is another source 
of irreversible (non-adiabatic) heating. 

As described in Reynolds et al. (2002) once the jet 
has ceased its activity and is turned "off" it loses its col- 
limation and subsequently collapses and disappears. The 
lobe of shocked jet material now rises into the heated at- 
mosphere due to its inertia and to its high entropy that 
makes it buoyantly unstable; the lobe evolves into a rising 
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Fig. 3. Time evolution of the M — 150 case. First row. density maps in logarithmic scale. Second row: specific entropy 
S = T jr?!'^ (logarithmic scale) which shows clearly the contribution of the entrainment by the high entropy material of 
the jet. Third row: entropy per unit volume (nlog(T/n^/'^)), which shows the contribution of the bow shock propagating 
into the medium. These quantities are shown in the columns at different times being r — 0.1 the time at which the jet 
is switched off and r = 3.9 the final simulated time. Notice the different spatial scale represented in the three columns. 
The density and temperature are given in units of their values in the center of the undisturbed medium. 



"plume" partly entraining the ambient medium contained 
in the shell. The bow shock, no more supported by the 
jet thrust, slows down and travels across the ambient as 
a shock wave or as a sound wave, depending on its en- 
ergy content compared with the energy of the unperturbed 
medium. 

These subsequent stages of evolution are represented 
in Fig. 121 where we show maps of some interesting quanti- 
ties at different times as they emerge from our simulation 
with M = 150: in the first row we represent density maps 
in logarithmic scale, in the second row we represent maps 



of entropy per particle {S = T/ri^/^) while, in the last row, 
maps of entropy per unit volume (5* = n log(T/n^/^)) are 
shown. The first column shows maps of these quantities 
at the end of the active phase (r — 0.1) while the other 
two columns show the evolution of the cocoon once the 
jet has ceased its activity until the final simulated time 
(t = 3.9). Notice the different spatial scale of the three 
columns. In the first row (density maps) we can visualize 
the different stages of evolution described above: during 
the active phase (left panel) the jets push against the sur- 
rounding medium that is compressed and heated in a thin 
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shell; in this phase the injected energy is far greater than 
the background energy and so the shock driven by the ex- 
panding lobe is very strong. Once the jets are turned off 
(central panel) the bow shock loses some of its strength 
but still expands as a mild shock wave since its energy 
is still greater than the background: the post-shock ma- 
terial is beginning to expand no more pushed by the ex- 
panding lobe. The lobes are stretching out due to their 
inertia and they are entraining the surrounding shocked 
ambient gas (even if the mixing effect was present yet in 
the active phase due to KH instabilities). In the advanced 
phases of evolution the lobes are evolving as rising plumes 
of low density and high entropy (see the second row) while 
the bow shock has lost quite completely its strength and 
so it mildly compresses the surrounding medium: in fact, 
at this spatial scale, the injected energy, that as we will 
see after is quite completely transferred to the ambient 
medium, becomes comparable to the background one (see 
the Ej/Eth ~ 0.4 ratio in Tabled). The shocked ambient 
material begins to recover a new hydrostatic equilibrium 
displacing the rising lobes as can be clearly seen in the 
rightmost panel of the third row too where the new core 
of higher entropy is clearly visible. 

In order to understand how much energy injected by 
the jet is transferred to the ambient and in which form, 
we have defined the following integral quantities 

Skin(t) - /r ^pv^ dV 

^pot(i) = jTp^dV ^^^^ 
Etotit) = Ei,,tit) + Ski„(t) + ^pot(0 

which represent respectively the internal, kinetic, gravita- 
tional potential and total energy of the ambient medium; 
in order to distinguish the ambient medium from the jet 
material, the integrals are performed weighing the ener- 
gies per unit volume with the tracer T defined as the sum 
of the four tracers % given in Sect. O and therefore mark- 
ing the entire cluster medium. The integrals are computed 
over the entire computational domain. Since the jet can 
transfer its energy to the ambient medium by means of 
the expansion work, done by the cocoon of thcrmalized 
jet material, and due to the mixing with its high entropy 
material at the contact discontinuity, it is also important 
to quantify the energy transferred to the ambient by these 
two processes separately. In the following discussion we 
will therefore indicate as "compressed" ambient medium 
the ambient material which forms the shell surrounding 
the cocoon: this medium has been compressed both adia- 
batically and irreversibly (at the shock front) by the ex- 
pansion work done by the cocoon. On the other hand, we 
will indicate as the "entrained" ambient medium the shell 
ambient material which is subsequently entrained by the 
cocoon. As noticed before in Reynolds et al. (2002) the 
ambient medium which is entrained by the cocoon can 
be distinguished from the compressed ambient medium 
since its specific entropy S = T/n'^/^, at least in the first 



phases of evolution, is much higher. Actually, the entropy 
per particle of the injected material is 10^ times higher 
than the entropy of the center of the unperturbed cluster 
medium, considering a jet in pressure equilibrium with the 
cluster core and with a density ratio v = 0.001 as in our 
main simulations. The entropy of the cocoon material is 
even higher since it is formed by the jet material which 
has passed through the terminal shock. On the other hand 
the bow shocks which are found in our simulations can not 
raise the entropy per particle of the cluster core to values 
greater than ~ 90 (in units of the central value) which is 
the value of the entropy per particle of the unperturbed 
cluster at r = i?vir. We therefore assume that the ambi- 
ent material contained in the central part of the cluster 
which has an entropy per particle S" > 5*0 = 90 has been 
entrained by the cocoon. This separation is not reliable at 
later times: when the bow shock reaches the outskirts of 
the cluster it raises the entropy to values S > Sq while nu- 
merical diffusion reduces the entropy of part of the cocoon 
below this threshold. We therefore perform this analysis 
only for times t/to < 0.4 taking 5*0 = 90 as entropy thresh- 
old to distinguish between compressed and entrained am- 
bient material. 

We can then define the total energy of the entrained 
(compressed) ambient material as: 

Etot. ontr (comp) = / T ( + -pt'^ + pcj)) dV (15) 

where E and C are the volumes in which S > Sq and 
S < So respectively. Clearly, Etot. cntr + -Btot, comp = Etot- 
We have plotted in Fig. 0] the fractional difference be- 
tween these ambient energies at time t and their initial 
values versus time. The values plotted are normalized with 
the value of the total injected energy. 




0.0 0.1 0.2 0.3 0.4 1 2 3 4 

t/t„ t/t„ 



Fig. 4. Temporal evolution of the variation of the am- 
bient energy. {Left panel) Temporal behavior of the to- 
tal ambient energy at the beginning of the cocoon evolu- 
tion {t/to < 0.4). We show plots of the total injected en- 
ergy (solid line), total ambient energy i?tot [dotted line), 
total energy acquired by the shocked compressed mate- 
rial i?tot, comp (dashed line) and by the entrained material 
-Etot. cntr (^dot-dashed line). (Right panel) Plotted here are 
the ambient internal energy Ei^t (solid line), the ambi- 
ent kinetic energy ii^kin (long-dashed line) and the ambi- 
ent gravitational potential energy i?pot (triple-dot-dashed 
line). Values of the energy are given in units of the total 
injected energy while time is given in units of Iq. 
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In the left panel of Fig. 21 we plotted the temporal 
behavior of the total injected energy (solid line) and the 
total energy transferred to the ambient medium i?tot (dot- 
ted line) in the initial stages of evolution (t/t^ < 0.4): it is 
possible to see that in the active phase the jet can transfer 
a high fraction of its energy to the surroundings 88% 
at t/to = 0.1). Plotted are also the increase of energy of 
the compressed ambient medium i?tot. comp (dashed line), 
which corresponds to the expansion work done by the co- 
coon, and the energy of the entrained ambient material 
Etot, cntr (dot-dashed line). As we can see at the end of the 
active phase (t/to — 0.1) the work done by the expanding 
cocoon corresponds approximatively to 60% of the total 
injected energy. This fraction agrees with the following 
simple analytical estimate (given also by Ruszkowski et 
al. 2003 in the case of an expanding bubble in pressure 
equilibrium with the surroundings). We assume the fol- 
lowing energy balance for the cocoon: 



dE,^^+pdV 
7-1 



(16) 



where £j is the energy injected in the cocoon, p is its av- 
erage pressure and V is its volume and so d(pV)/('y — 1) 
is the variation of its internal energy and pdV is the ex- 
pansion work done at the contact discontinuity. Assuming 
also the spherical self-similar behavior given by Eqs. (|ll|l 
and H12|) for the injected energy, the atmosphere profile 
and the cocoon radius, we obtain the following efficiency e 
(ratio between the expansion work and injected energy): 



pdV _ (7-l)(3d + 6) 
dEj" ~ ^(3d + 6) + dk + 2d-6 



(17) 



which for d — 1 (constant injection power) and k — —2 
(isothermal sphere profile) gives e = 66%. 

At t/to = 0-1 the energy transferred to the ambient 
medium by means of the mixing with the cocoon high en- 
tropy material corresponds then to the remaining 28% 
of the injected energy but it is possible to see that, just 
after the jet stopped its activity, the entrained material 
tends to decrease its (mostly thermal) energy due to the 
expansion of the cocoon which is no more fueled by the 
jet. On the other hand, the ambient medium compressed 
by the still expanding cocoon increases its energy up to 
values ~ 90%i?j. As pointed out before, it is difficult to 
distinguish between the compressed and the entrained am- 
bient medium on longer timescales, but based on these 
considerations, we can argue that at later times most of 
the injected energy is retained by the compressed ambient 
medium and only a small fraction by the rising lobes and 
by the entrained ambient material. 

Looking at the right panel of Fig. 0] we can finally see 
in which form the energy is transferred to the ambient 
medium on longer timescales. In this plot we make no dis- 
tinction between compressed and entrained material. In 
the active phase of the jet (t/to < 0.1) the ambient is com- 
pressed, heated by the expanding shock wave and partly 
entrained by the cocoon, thus increasing its internal E'int 



t/t. 




Fig. 5. (Left panel) Temporal evolution of the average 
entropy per particle S—T/n'^^^ of the ambient material 
marked by the innermost three tracers: 7i (r/vg G [0, 1], 
dashed line), T2 (r/r^ £ ]lj3], dot- dashed line) and Ta 
(r/vs G ]3, 6], solid line). (Right panel) Temporal evolution 
of the total entropy increase Stot (Eq. (|18|l l of the inner- 
most tracer 7i: plotted are the total entropy increase of 
the ambient medium (dashed line), the entropy increase 
of the shocked material (triple-dot-dashed line) and the 
entropy increase of the entrained material (dotted line). 
Time is given in units of t^. 



(solid line), kinetic i^kin (long-dashed line) and gravita- 
tional potential energy E'pot (triple-dot-dashed line) since 
its material is pushed higher in the potential well. In the 
passive phase the ambient material, no more pushed by 
the internal cocoon, slows down, decreasing its kinetic en- 
ergy. Even if the bow shock is still expanding and heating 
the ambient, behind it the shocked material can expand 
and decrease its internal energy trying to recover the hy- 
drostatic equilibrium. Only the gravitational energy is still 
growing since the shock still pushes the material higher in 
the potential well while the expanding shocked material 
tries to recover the gravitational equilibrium on a higher 
adiabat thus increasing its potential energy (it is less dense 
than it was in the initial equilibrium). 

Since not all the energy that is transferred to the intr- 
acluster medium is dissipated, what it is more important 
to determine is the energy that is irreversibly transferred 
into the ambient medium or rather the heating that de- 
termines the increase of the entropy of the ICM. As we 
have pointed out in this section the entropy of the am- 
bient medium can be raised in two ways: shocked by the 
blastwave pushed by the expanding cocoon or through the 
interaction with the high entropy material of the rising 
lobes. In the following Section we will quantify the energy 
transfer associated with these two phenomena. 

5. Entropy evolution 

We can understand how the ambient entropy has increased 
due to the interaction with the jets looking at the second 
and third row of Fig. 13 the maps in the second row (evolu- 
tion of the entropy per particle S — T jr?!"^) show clearly 
how the entrainment with the high entropy material of 
the jets can raise substantially the entropy of the mixed 
material while the entropy increase of the shocked gas is 
not even noticeable. But if we look at the entropy per unit 
volume (third row) we see clearly the expanding shock to 
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be the dominant feature: this means that the entrainment 
with the high entropy material of the lobes increases the 
entropy of the ambient to higher values but involves only 
a little fraction of the mass of the ICM; on the other hand 
the expanding shock wave increases the ambient entropy 
to lower values involving a greater fraction of the ICM 
traveling across a big portion of the structure. This be- 
havior can be seen more quantitatively in the left panel 
of Fig. where we have plotted the evolution of the mean 
entropy per particle of the first three tracers that follow 
the evolution of the ambient medium: we can see that the 
entropy increase of the tracers is mostly associated with 
the passage of the bow shock since once the shock has 
overtook the inner tracer this one evolves adiabatically 
while the outer one starts to increase its entropy due to 
the shock passage. To evaluate the importance of entrain- 
ment in rising the entropy of the ambient medium we can 
define separately the variation of the total entropy of the 
entrained material, characterized by an entropy per par- 
ticle S > Sq — 90, and of the compressed material, which 
has S < Sq. We remember that this kind of analysis is not 
correct on longer timescales. We plotted in the right panel 
of Fig. 121 the variation of the total entropy 



S'tot ^ Tin log 



T 

72J3 



dV 



(18) 



of the inner tracer at the beginning of the cocoon evolution 
{t/tf) < 0.4) (dashed line, values are normalized to the 
plateau reached by this curve). Computing the integral 
of Eq. (|18|l over the volumes in which Tjr?!'^ < Sq and 
T jr?!'^ > 5*0, we plotted respectively also the variation of 
entropy of the ambient material compressed by the bow- 
shock [trvple- dot- dashed line) and the variation of entropy 
of the entrained material [dotted line): we can see that the 
shocked material can account for a high fraction of the 
total increase of entropy (~ 94%) while the contribution 
of the entrained medium is negligible (~ 6%). 

On this basis we can argue that in our simulations the 
shock heating is the dominant effect determining the ther- 
modynamical evolution of the ICM. Since the propagating 
shock, mostly after the jet has been turned off, assumes 
a spherical shape, we can define an entropy profile taking 
averages over spherical shells of given mass, building an 
entropy profile as a function of the mass enclosed inside a 
given radius. In Fig.|Slwe plot the entropy profile associ- 
ated with the passage of the shock wave as a function of 
the enclosed mass for the six cases of Table ^ the three 
panels represent the three clusters taken into account (Tow 
= 2, 1, 0.5 keV from left to right) while the solid lines in 
each panel represent the entropy profiles determined by 
the two jet powers considered (10^^ and 5 x 10^^ erg s~'^ 
for the lower and upper curves respectively). The unper- 
turbed initial entropy profiles are shown in each panel with 
a dotted line. 

In the cases in which the ambient energy enclosed in- 
side the shock radius becomes comparable to the almost 
constant energy associated with the shock wave, the shock 
loses its strength and evolves slowly as a sound wave, 
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Fig. 6. Entropy S = TjrP'^'^ profiles as a function of 
Lagrangian (mass) coordinate at the beginning and at the 
end of the simulation. The panels refer to the simulations 
performed with the Tew — 2 keV (lejt •panel), Tew = 1 keV 
(central panel) and Tow = 0.5 keV (right panel) clusters. 
In each panel are shown the initial entropy profile (dotted 
line), the profiles determined by the interaction with the 
jets characterized by a luminosity Lj = 10^^ erg s~^ (lower 
solid line) and Tj = 5 x 10^^ erg s~^ (upper solid line) and 
analytical extrapolations of the profile up to the virial ra- 
dius (dashed line). The mass coordinate is given in units 
of the total gas mass content of the clusters /bar-^vir- 

not increasing the ambient entropy from a certain radius 
(mass) onward. From this point onward the entropy pro- 
file nearly matches that of the unperturbed medium (see 
for example the lower solid lines in the left and middle 
panel in Fig. 1^. For the cocoons that remain overpres- 
sured for the entire evolution inside the cluster or group, 
the expanding shock never loses its strength and the en- 
tropy of the atmosphere is increased everywhere. In some 
of our simulated cases, we cannot follow the evolution of 
the blastwave up to the the virial radius, since the virial 
radius is not contained inside the computational domain 
of some simulations, or it is not reached at the final time 
of the simulation. In Fig.|Slthe entropy profiles derived di- 
rectly from the numerical simulations (solid lines) stop at 
the farthest coordinate values reached. In this situations, 
we build an extrapolation of the entropy profile from the 
maximum radius reached in the simulation up to the virial 
radius in the following way: we assume that in these final 
stages of evolution the energy E associated with the blast 
wave is a constant; we assume that a constant fraction of 
this energy goes into thermal energy (pressure) and that 
this pressure (> external pressure) drives the expansion 
of the shock. We can then write the following expressions 
for the pressure Ps driving the shock and the radius rg of 
the shock (see Appendix A in Zanni et al. 2003): 



E 



dt 



P{rs) 



where p(rs) is the ambient density evaluated at the shock 
radius rg. This set of equations gives the following solution: 



t-tn^K 



p(r) ^r^dr 



(19) 



where tq and to are the shock radius and the time at which 
we begin to extrapolate the entropy profile while if is a 
constant that is determined matching the solution for the 
averaged shock radius coming from the simulations with 
the solution of Eq. ifTi^ . Once we have determined the 
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shock radius as a function of time we can derive its speed 
and its Mach number with respect to the ambient sound 
speed. From the Mach number of the shock we can finally 
derive the entropy jump associated with the expanding 
wave. These extrapolated profiles are represented in Fig. 
Elwith dashed lines. 

From these entropy profiles it is possible to define an 
"equivalent heating" assuming that these profiles are ob- 
tained from an isochoric transformation, i.e. keeping con- 
stant the initial density profile of the ambient and increas- 
ing only its temperature. This quantity gives an informa- 
tion on the amount of energy that is irreversibly dissipated 
inside the cluster by the jet and goes into heat. This es- 
timate is obviously different from the real amount of in- 
jected energy that is converted into heat, since the ther- 
modynamical transformation associated with the propa- 
gation of the shock which is compressing the gas is not 
isochoric. Anyway this is exactly the amount of energy 
gained by the gas particles in a non-adiabatic compression 
followed by an adiabatic expansion to recover the initial 
density. In Fig. [3 we show a radial profile (in mass coor- 
dinates) of this equivalent heating, showing the different 
amount of heating energy per particle released inside 
the clusters for our six cases: in the three panels we show 
the results for the three clusters considered (Tew = 2, 1, 
0.5 keV from left to right), and in each panel the solid 
lines represent the profiles obtained with the two jet pow- 
ers taken into account (10"*^ and 5 x 10"*^ erg s~^ for the 
lower and upper curves respectively ). The dashed lines 
represent the heating derived from the extrapolated en- 
tropy profiles. In Fig.|Slwe show the temporal behavior of 
the spatially integrated value of this quantity normalized 
to the total injected energy, showing that in all the six 
cases considered a high fraction (up to ^ 75%) of the in- 
jected energy is dissipated through irreversible processes 
and heats the cluster. In Fig. 01we showed that almost all 
the injected energy is transferred to the ambient medium: 
we find a lower "equivalent heating" since it accounts only 
for processes which can raise the entropy of the ICM and 
therefore neglects the energy which is transferred adiabat- 
ically. The upper curves in Fig. |H1 refer to the cases with 
a E-jEt\i ratio equal to 10.9, 1.9, 2.2, the middle curves 
to the cases characterized by Ej/Eth = 0.39, 0.34 while 
the lower solid line represents the £j/i?th — 0.068 case. 
The jets characterized by a low Ej/Eth ratio can also ef- 
ficiently dissipate their energy, since they can be in an 
overpressured phase at the beginning of their evolution. 
Observations of the diffuse X-ray thermal emission asso- 
ciated with the low power radio source NGC1052 (Kadler 
et al. 2004) could confirm this possibility. 

6. Hydrostatic equilibrium 

Once we have calculated the entropy profiles as a function 
of the included mass, assuming that the atmosphere will 
recover the hydrostatic equilibrium without additional dis- 
sipation, we can determine the radial behavior of all the 
thermodynamical quantities in the new hydrostatic equi- 
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Fig. 7. Plot of the heating energy per particle i?h needed 
to obtain the entropy profiles of Fig. as a function of 
mass coordinate. The panels refer to the Tow = 2 keV {left 
panel), Tow — 1 keV {central panel) and Tow = 0.5 keV 
{right panel) clusters. In each panel we plot the heating 
energy per particle given by the Tj = 10"*^ erg s~^ {lower 
solid line) and ij = 5 x 10^^ erg s^^ jets {upper solid 
line) and analytical approximations up to the virial ra- 
dius {dashed line) obtained from the extrapolated entropy 
profiles. The mass coordinate is given in units of the total 
gas mass content of the clusters /barAfvir- 
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Fig. 8. Plot of the total heating energy dissipated by the 
jet in units of the total injected energy as a function 
of time. The curves refer to the cases characterized by 
M = 105 {solid line), M = 150 {dotted line), M = 211 
{dashed line), M =180 {dot- dashed line), M = 255 {dou- 
ble dot-dashed line) and M = 361 {long dashed line). It 
is important to notice that the lines are divided in three 
groups characterized by different ratios Ej/Eth between 
the total injected energy and the thermal energy of the 
cluster: Ej/Eth = 10.9, 2.2, 1.9 (upper lines), Ej/Eth = 
0.39, 0.34 (middle lines) and Ej/Eth = 0.068 (lower line). 
See Table ^ for the correspondence between Ej /Eth and 
Mach number. 

librium solving the following system of equations: 
dP - ( P \^ „ 

dr — \S(M) J y 



< 




where we used the relation S{M) = P/p'* , where S{M) is 
the averaged entropy profile as a function of the included 
mass as calculated from our simulations. The gravity ac- 
celeration g is derived from the potential of Eq. ((SJ. In 
order to solve this system of equations we need two bound- 
ary conditions being the first, obviously, M{r = 0) = 0. 
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The second boundary condition is given requiring that 
the thermal pressure found in correspondence of the virial 
(gas) mass is the same given in the initial profile, assum- 
ing, in a sense, that the conditions outside the cluster 
did not change. However, we found that the quantities 
that we calculated from these final profiles did not depend 
strongly on the assumed boundary conditions: our outer 
boundary condition requires only that the pressure outside 
the cluster is much lower than the central one. The cen- 
tral values of the interesting thermodynamical quantities 
(density, pressure and temperature) depend only slightly 
on the boundary condition assumed, while the integral 
quantities such as the gas mass contained inside the virial 
radius, as noticed before by Voit et al. (2002), is more sen- 
sible. As an example we plot in Fig. |51the density profiles 
of the new hydrostatic equilibrium for the six simulations 
of Tabled in the three panels we show the results for the 
three clusters considered (Tew = 2, 1, 0.5 keV from left to 
right), the dot-dashed and the dashed lines represent the 
profiles obtained with the 10^^ and 5 x 10"*^ erg s~^ jets 
respectively. 
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Fig. 9. Radial density profiles of the clusters defined in 
the initial conditions (Eq. (jSJ) and obtained by calculat- 
ing the hydrostatic equilibrium determined by the entropy 
profiles of Fig. El The panels refer to the Tow = 2 keV 
{left panel), Tow = 1 keV (central panel) and Tow = 0.5 
keV (right panel) clusters. Plotted are the profiles of our 
baseline model (solid line) and the profiles determined by 
taking into account the energy dissipated by the Tj = lO''^ 
erg s~-^ (dot-dashed line) and Lj = 5 x 10^^ erg s~^ jets 
(dashed line). The density is given in units of the central 
density in the initial conditions while the radial coordinate 
is normalized with the characteristic radius r,. 



Another way to give an estimate of the energy per par- 
ticle that the jet has dissipated into the ambient medium is 
to evaluate the difference between the average energy per 
particle (thermal and potential) in the final and the ini- 
tial hydrostatic conditions. Moreover, this estimate does 
not depend on the assumption of an isochoric transforma- 
tion on which the calculation of the "equivalent heating" 
is based. In Fig. ^] we plot in the left panel the excess en- 
ergy Ecxc of the particles contained inside the virial radius 
in the final state versus the temperature of the cluster due 
to the 10*^ erg s^^ jets: the total energy is shown with the 
crosses while the asterisks and the diamonds correspond 
respectively to the potential and thermal excess energy. 
We can see that in the final state the greater amount of 
excess energy is in the form of potential energy since the 
heating of the ICM has led to a huge expansion of the 



medium thus increasing its gravitational energy (the final 
state is less dense then the initial one and the gravitational 
potential is the same). The expansion of the structure has 
obviously led to a substantial adiabatic cooling and so the 
final thermal energy is not much greater than the initial 
one. In the right panel of Fig. ^| we plot the total ex- 
cess energy as a function of the mean temperature of the 
clusters as due to the 10"*^ (crosses) and 5 x 10"*^ erg s~^ 
(asterisks) jets. We can see that the excess energy is lower 
than the energy per particle dissipated in the computa- 
tional domain (Table |3J|: this is due to the work done in 
the expansion to move the mass that now is located out- 
side the virial radius and so this difference is more evident 
in the cases that have expanded more. 




Fig. 10. Plot of the difference between the average energy 
per particle of the final and initial hydrostatic equilibrium 
as a function of the final emission- weighted temperature of 
the clusters. The average of the energy excess is calculated 
over the particles contained inside the virial radius in the 
final hydrostatic conditions. (Left panel) Plot of the aver- 
age energy excess per particle determined by the Lj = 10^^ 
erg s~^ jets: plotted are the total energy (crosses), the 
gravitational potential energy (asterisks) and the internal 
energy (diamonds). (Right panel) Plot of the total energy 
excess per particle determined by the Lj = 10"*^ erg s~^ 
(crosses) and the Tj = 5 x 10^^ erg s^^ jets (asterisks). 



With these new profiles it is possible to determine how 
the observable properties of the clusters have changed 
consequently to the jet energy injection. We then com- 
puted the X-ray luminosity adding to the bremsstrahlung 
emissivity a correction for metal-line cooling assuming a 
plasma with Z — 0.3 Zq metallicity. In Fig. ^2 we plotted 
the results obtained for the X-ray luminosity - tempera- 
ture relation while the comparison between the observed 
and the simulated Lx — M200 relation is shown in Fig. 

In both Figures we also plotted the scaling of our ini- 
tial model for the X-ray luminosity including line emis- 
sion (solid line) and for bremsstrahlung luminosity only 
(dashed line): we recall that this last quantity satisfies the 
self-similar behaviors Lbr oc To^w ^-^i^ Tbr oc Mggg respec- 
tively. It is possible to see that the luminosities obtained 
with the higher power Lj = 5 x 10''^ erg s~^ jets are in 
good agreement with the observational data, mostly at 
the smaller simulated scales. From Fig. we can also no- 
tice that, as discussed before, the temperature of the final 
state is only slightly higher than the initial one and the 
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Fig. 11. Bolometric luminosity-emission weighted tem- 
perature relation for our simulations characterized by 
Lj = 10^^ erg s~^ {empty squares) and Lj = 5 x 10''^ 
erg s~^ (crosses) in comparison with observational data 
(pluses: Helsdon & Ponman 2000; filled squares: Arnaud & 
Evrard 1999; filled triangles Markevitch 1999). Plotted are 
also the scaling of the bremsstrahlung luminosity (dashed 
line) and the X-ray luminosity with a correction for line 
emissivity (metallicity 0.8^0) (solid line) of our initial 
model. The bolometric luminosities of the three initial un- 
perturbed clusters (asterisks) and of the four Tornatore et 
al. (2003) simulations (diamonds) from which our initial 
conditions were derived are shown for comparison. 
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Fig. 12. Comparison between the X-ray luminosity - 
-^^200 relation obtained from our simulations with Lj — 
10'^^ erg s~^ (empty squares) and Lj = 5 x 10*^ erg s""'^ 
(crosses) and the observational data (triangles: Rciprich 
& Bohringcr 2002). Plotted are also the scaling of the 
bremsstrahlung luminosity (dashed line) and the X-ray lu- 
minosity with a correction for line emissivity (metallicity 
O.SZq) (solid line) of our initial model. The bolometric 
luminosities of the three initial unperturbed clusters (as- 
terisks) are also plotted. 



the observed M200 data represent the gravitational mass 
contained inside -R200, we determined M200 in our simula- 
tions adding up the dark matter content and the gas mass, 
which is not self-gravitating. Due to the low gas mass frac- 
tion of these structures this should not affect dramatically 
our results. Finally, since neither the average temperature 
nor the total mass content of the clusters are strongly af- 
fected by the energy injection, the Tcvj—M2oo relation does 

2/3 

not deviate much from the Tow oc M2QQ self-similar behav- 
ior, consistently with the observations (see Horner et al. 
1999 and Reiprich & Bohringcr 2002). In Fig.HSlwe plot 
the relation between the entropy per particle estimated at 
0.1i?20O and the mean temperature. As previously noticed 
for the X-ray luminosity relations, also the increase of en- 
tropy is affected more by the lower density of the final 
state than by the slightly higher temperature. In this plot 
is even more evident that a good agreement with the data 
is obtained only with the higher power energy injections 
at the smaller scales (Tew < 1 keV) . This in not surprising: 
it has been argued by Roychowdhury et al. (2004) that an 
energy injection which increases linearly with the cluster 
mass is needed to reproduce the entropy profiles observed 
by Ponman et al. (2003): on the other hand we are inject- 
ing the same amount of energy at all scales which does 
not seem to be enough for the more massive clusters. 
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Fig. 13. Entropy (at 0.1i?2oo) - emission weighted tem- 
perature relation for our simulations (Lj ~ W^^ erg s~^, 
empty squares and Tj = 5 x 10"*^ erg s^^, crosses), in com- 
parison with observational data (pluses: Ponman et al. 
2003). The observational data represent averages in dif- 
ferent bins of temperature. The dotted line represents the 
scaling of the initial model (Eq. (jSJ). The dotted-dashed 
line represents the constant floor 125 (ft./0.7)^^/^ keV cm^ 
(crosses) determined by Lloyd-Davies et al. (2000). 



lower X-ray luminosity is therefore mainly determined by 
the lower density of the final state. 

On the other hand (see Fig. 1121) . Af200 is poorly af- 
fected by the gas mass loss, since the mass content of the 
cluster is mainly determined by the dark matter. While 



7. Tests on numerical resolution and jet 
parameters. 

We finally discuss the results obtained with the "test" 
cases reported in Table |21 We recall that these simula- 
tions have been performed in order to test how much the 
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results obtained with our six simulations depend on the 
grid resolution and on the jet parameters assumed; we 
therefore took a reference case with an initial atmosphere 
with Tew = 0.5 keV (M = 211, third hne in Table [H and 
perform three additional runs. The first one (first line in 
Table |5} is a case with the same jet parameters (Mach 
number, jet density, jet radius and duration of injection) 
of the reference simulation but with half the grid resolu- 
tion: we don't find significant differences in the entropy 
profiles (and therefore in derived results), showing that 
the adopted resolution is high enough to obtain numeri- 
cal convergence. On the other hand, it is well known that 
shock capturing methods such as the PPM rapidly con- 
verge to a solution in problems involving strong shocks 
as in our simulations. The case represented in the second 
line of Table|21has the same jet kinetic power and injection 
time of the reference simulation but higher Mach number 
and lower density: also in this case, discrepancies between 
the resulting entropy profiles and the ones described in 
the previous Sections are well within a few percent. This 
indicates that the parameter controlling the jet dynamics 
is the jet power, quite independently from other parame- 
ters (Mach number, density and radius) which character- 
ize the jet. This hypothesis is confirmed also by the fact 
that simulations characterized by similar ratios between 
the injected energy and the thermal energy of the unper- 
turbed cluster Ej/Eth but with different Mach numbers 
and different ratios between the jet radius and the core 
radius r-Jr^, have given similar results: it has been yet no- 
ticed that similar efficiency curves (Fig. IHl are obtained 
from simulations with comparable Ej/Eth- Moreover Fig. 
El shows clearly that these simulations give similar final 
density profiles with approximatively the same decrease in 
the central density (compare for example the dashed line 
in the central panel with the dot-dashed line in the right 
panel which are obtained from two simulations character- 
ized by Ej/Eth — 1.9 and 2.2 respectively): this means 
that if these cases are rescaled to the same initial con- 
ditions (same Eth) they will give comparable observable 
quantities. 

Different considerations can be done for the last case 
in Table El (third line) that is a jet with half the kinetic 
power of the reference simulation (due to its lower den- 
sity) but injected for twice the time {t/to = 0.2) so as to 
have the same total injected energy. Since the strength of 
the expanding shock is mainly determined by the injected 
energy, it is easy to guess that in the central part of the 
group the shocks driven by the lower power case will be 
weaker, while on a larger scale, since the energy injected 
by the two jets is the same, the shocks will have more or 
less the same strength. This is confirmed by Fig. 1141 where 
we plot the entropy profiles given by the lower power (solid 
line) and the higher power jet (dashed line) as a function 
of the mass coordinate: while toward the central part the 
lower power jet determines a lower entropy profile, on a 
larger scale the two profiles are absolutely similar. Also 
on the large scale the little discrepancy that can be seen 
between the two profiles can be ascribed to the small dif- 
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Fig. 14. Entropy S = T jn}^'^ profiles as a function of 
Lagrangian (mass) coordinate. Plotted are the initial con- 
ditions (dotted line) and the results obtained with the 
M = 211 simulation (dashed line) which represents a 
ij = 10**^ erg s~^ jet injected in a Tew = 0.5 keV clus- 
ter for a time t ^ 2 x 10^ years and with the case in the 
third line of Tahle^ (solid Zme)which is a jet with half the 
kinetic power (lower density) injected for twice the time. 
The mass coordinate is given in units of the total gas mass 
content of the clusters /bai-Afvir- 



ference between the two total injected energies. The two 
jets inject the same amount of kinetic energy but, having 
the same internal energy flux, the lower kinetic power jet 
injects twice the internal energy, giving a difference of a 
few percent on the total injected energy. The difference 
found in the central part of the profiles implies therefore a 
lower (^ 30%) "equivalent heating" and leads obviously to 
a different new equilibrium, with a higher central density 
and lower temperature. This produces a higher (^ 30%) 
X-Ray luminosity. 



8. Discussion and conclusions 

In the previous Sections we have shown how a (double- 
sided) powerful jet can dissipate efficiently its bulk kinetic 
energy thus rising the entropy of the medium through 
which it propagates. This was done by analyzing by means 
of numerical hydrodynamical simulations the interaction 
of two powerful supersonic jets with different powers (Lj = 
10^^, 5 X 10**^ erg s~^) with the ICM of three gravitation- 
ally heated clusters of galaxies characterized by different 
temperatures (Tqw = 0.5, 1, 2 keV). 

The mechanism through which these jets interact with 
the ambient medium can be summarized as follows: an 
underdense supersonic jet is slowed down by a terminal 
shock (Mach disk) where its kinetic energy is thermalized 
and inflates an overpressured cocoon; the expansion of the 
cocoon compresses the surrounding medium and, as long 
as the energy transferred to the ambient by the expansion 
work is higher than the background energy (even after the 
jet has ceased its activity), the compression of the medium 
will be mediated by an expanding shock (bow shock) (see 
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Eq. (|13|l ') where the energy transferred to the medium can 
be dissipated. 

A high fraction of the injected energy can be deposited 
into the ambient medium by the expansion work (~ 60% 
at the end of the active phase but even higher afterward, 
see Fig. 0} while only a small fraction is retained in the 
cocoon and in the entrained ambient material (< 10%). 
As it has been shown e.g. by Ruszkowski et al. (2003), an 
expanding bubble in pressure equilibrium with the sur- 
roundings can transfer to the medium around 40% of the 
injected energy but, since mostly in the active phase of in- 
jection the expanding cocoons are strongly overpressured 
this fraction can be much higher in our simulations. 

Once the jet has ceased its activity it quickly loses 
its coUimation and disappears while the high entropy co- 
coon formed by the thermalized material of the jets buoy- 
antly rises into the atmosphere entraining the ambient 
medium. As proposed in other works the energy of the 
rising plumes (or bubbles) associated with the relics of 
the radio lobes can be dissipated thanks to the expansion 
work done as they rise in the decreasing pressure ambient 
(Roychowdhury et al. 2004), to the mixing between high 
and low entropy material (e.g. Dalla Vecchia et al. 2004) 
or to the viscous dissipation of sound waves generated by 
the rising lobes (Ruszkowski et al. 2004) and generally is 
high enough to balance cooling losses and limit cooling 
flows in the cluster cores. But since we have shown that 
the energy retained in the rising plumes is negligible com- 
pared to the energy dissipated at the shock front, in our 
discussion we neglected the effects due to the entrainnient 
between the high entropy material of the plumes and the 
ambient medium. 

The energy dissipated at the shock front has been eval- 
uated looking at the entropy profiles of the shocked ambi- 
ent medium: we found that after the passage of the bow 
shock the entropy profile as a function of Lagrangian mass 
coordinates tends to reach a steady condition and only 
small fluctuations due to the action of the rising lobes are 
observed, thus confirming that in our simulations the ex- 
panding bow shock is the main heating mechanism. As it 
was shown in Fig. El the entropy of the ICM can be sub- 
stantially raised as long as the energy associated with the 
expanding shock is greater than the energy of the unper- 
turbed medium contained inside the shock radius: as can 
be seen in Fig. |B1 (upper solid curves) and Fig.ll3lonlv the 
higher power jets (5 x lO''^ erg s~^) can raise the entropy of 
the central core of clusters up to values > 125 (/i/0.7)^^/^ 
keV cm^ which represents the "entropy floor " determined 
by Lloyd-Davies et al. (2000). Moreover, the data obtained 
with the higher power jets are in good agreement with the 
average core entropies determined by Ponman et al. (2003) 
only for the small groups (Tew < 1 keV). 

A huge amount of energy is required in order to obtain 
the observed core entropy with a single energetic event. 
Using the same isochoric approximation adopted to cal- 
culate the "equivalent heating" , we estimate that average 
energies around 5.3 keV, 4.6 keV and 3 keV per particle 
for the Tew = 0.5 keV, Tew = 1 keV and Tew = 2 cluster 



Table 3. Average energy per particle dissipated by jets of 
different power Lj in the three clusters characterized by 
Tew=0.5, 1, 2 kev. Reported are the average energy per 
particle dissipated inside the entire domain (first value) 
and inside the cluster core R < Vs (second value) at the 
end of the simulations. 
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respectively are required to increase the average core en- 
tropy of our baseline model up to values > 125 (/i/0.7)~^/^ 
keV cm^, with a single energy injection at z = 0. In Fig. 
Owe plotted the equivalent amount of energy per particle 
required to increase the entropy of the initial atmosphere 
up to the entropy determined by the jet injection while in 
Table 121 we show for every simulation the average energy 
per particle dissipated by the jets inside the whole com- 
putational domain (first value) and inside the cluster core 
(second value) at the end of our simulations: these values 
are calculated averaging the values of the curves plotted 
in Fig. [71 over the mass contained in the whole compu- 
tational domain and in the cluster core respectively. We 
can see that the energy dissipated by the 5 x lO'*^ erg s~^ 
jets in the three cluster cores is high enough to raise the 
core entropy up to the entropy floor but nevertheless the 
lower power cases can give at least in the cluster cores the 
heating energy around 1 keV per particle required in cos- 
mological simulations to satisfy observational constraints. 
The low power cases can efficiently dissipate their energy 
inside the cluster core since at least at the beginning of the 
evolution their cocoons are overpressured with respect to 
the ambient medium, thus driving strong shocks through 
it. On a larger scale these shocks loose their strength, 
and they thus resemble observations of sound waves and 
weak shocks propagating in the ICM (Fabian et al. 2003, 
Forman et al. 2003). 

The duration of the overpressured phase during which 
the jets can efficiently dissipate their energy is mainly de- 
termined by the ratio between the injected energy and the 
energy of the unperturbed medium (sec Table ^ and so 
it depends strongly on the characteristics of the medium 
through which the jets propagate. As shown by Reynolds 
et al. (2002) even if a jet can efficiently shock the core of 
a rather shallow entropy profile (they assumed an isother- 
mal /3-model with /3 = 0.5) the lower entropy material out- 
side the shocked region will flow inward taking the place 
of the strongly shocked medium thus minimizing the effect 
of the energy dissipation. On the other hand the steeper 
entropy profiles that we assumed as initial condition as 
representative of clusters whose ICM is only heated by 
gravitational processes are affected more significantly by 
the interaction with a supersonic jet. 
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The fraction of the injected energy that is dissipated 
inside the cluster is plotted in Fig. |Slfor every simulated 
case as a function of time. In every case, a high fraction (up 
to 75%) can be irreversibly dissipated but a few differences 
can be noticed: for instance, a rather small dependency 
of this efficiency on the ratio Ej/Eth between the total 
injected energy and the thermal energy of the cluster (with 
« 10% differences). Cocoons characterized by a higher 
ratio can remain in an overpressured phase for a longer 
time driving shocks across a higher fraction of the ICM 
and dissipating more efficiently their energy. As noticed 
before, also less powerful jets can dissipate their energy 
onto the ICM, at least at the beginning of their evolution. 
Therefore, many less powerful jets should be able to heat 
effectively the ICM as much as one single powerful event. 

Taking into account the simulated case characterized 
by Ej/Etu = 10.9 (Lj = 5 X lO-^^ gj.g s-\ Tew = 0.5 keV) 
which is one of the cases that better satisfies the observa- 
tional constraints (see Figs. II II and 113(1 . it is evident that 
at the Tew = 0.5 keV scale an energy around 8.7 times 
greater than the initial thermal energy of the unperturbed 
cluster (assuming an efficiency around 80%) must be dis- 
sipated in order to heat the ICM up to the observed val- 
ues. This estimate is consistent with the results obtained 
by Roychowdhury et al. (2004) even if they modeled the 
AGN activity in a different way (they considered the work 
done by a steady flux of buoyant bubbles) : their linear re- 
lation between the cluster mass and the energy that must 
be dissipated in the ICM in order to match the entropy 
excess at 0.1i?20O and -R500 observed by Ponman et al. 
(2003) is consistent with our results on the overall energy 
budget at the Tow = 0.5 keV scale. A somewhat different 
estimate was given by Cavaliere, Lapi & Menci (2002): to 
match the observations of the L — T relation for a T = 0.5 
keV cluster they require an energy which is around 3 times 
the energy of the unperturbed medium. They modeled the 
AGN activity with a self-similar mild blastwave propagat- 
ing through the ICM pushed by a piston which injects 
energy in the ICM with low power but for long dynamical 
timescales: such a blastwave does not efficiently heat the 
surrounding medium but tends to push the ICM out of 
the virial radius thus decreasing the mass content of the 
cluster and then its X-ray luminosity. In our simulations, 
clusters are heated more efficiently, and the impact of this 
on the L — T relation can be evaluated. Heated clusters 
tend to expand and cool down adiabatically: when they re- 
cover the gravitational equilibrium they have a less dense 
core and a temperature only slightly higher than the ini- 
tial one thus decreasing their X-ray luminosity. As it is 
shown in Fig. ^] the cases that are heated up to the ob- 
served values (Fig. 113(1 can expand to a new equilibrium 
to give X-ray luminosities which are consistent with the 
observed L — T relation. 

From the above discussion, our conclusion is that the 
mechanism discussed here can represent one viable way to 
induce a non-gravitational heating of the ICM. An extra 
heating of « 1 keV per particle, is known to be sufficient 
for bringing the theoretical Lx — T relation in agreement 



with observations (see e.g. Borgani et al. 2002). We have 
shown how such a heating level can be provided by AGN 
jets having kinetic power in the range lO''^ — 5 x lO''^ erg 
s-i. 

The situation of one single, powerful jet which heats 
the ICM to the observed level can be quite rare in real- 
ity. More probably, several, less powerful events contribute 
to the heating. We showed that the contribution of such 
events, while smaller, have a non negligible effect. Other 
non gravitational mechanism, like the "effervescent heat- 
ing" or other interactions between AGN-originated bub- 
bles or plumes of gas and the ICM, could also be simulta- 
neously in action, and we regard them as alternative, and 
not exclusive, respect to the one described here. 

Another important detail is that the calculation re- 
ported here are based on the ICM physical characteristics 
at a redshift z = 0, while most probably a jet heating 
would occur at higher redshifts, in the environment given 
by proto-cluster structures. One could speculate that, de- 
pending on the actual number and power of heating events 
which happened in these proto-clusters, the final excess 
heating can significantly vary among different clusters, 
thus explaining not only the break in the self-similarity, 
but also the observed spread in the X-ray properties at a 
given temperature. For instance, from Fig. II II one can see 
that, at a temperature T ^ 0.7 — 0.9 keV, the observed 
luminosity can vary by more than an order of magnitude. 
If the extra-heating of the clusters atmospheres is due to 
the AGN activity, and if such activity is somehow linked 
with the proto-cluster "cosmological" environment, the 
spread could be reminiscent of the formation history of 
the cluster. 

To put on firmer ground such speculations, we plan 
to repeat the current analysis using several atmospheres 
of proto-cluster structures taken from N-Body/SPH sim- 
ulations at high redshifts. This will allow us to quanti- 
tatively asses the contribution of AGN jets heating on 
the gas which will form the clusters atmospheres at red- 
shift z = 0, and possibly to follow the heated gas during 
the cluster formation process, again with the use of N- 
Body/SPH simulations. This phase is preliminary to the 
attempt to give a parametrization of the effect of AGN 
activity on the ICM, or at least of the contribution due 
to the jets heating, which could be directly used in self- 
consistent cosmological simulations. 
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